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We have investigated the behavior of the resistance of graphene at the n — Landau Level in 
an intense magnetic field H. Employing a low-dissipation technique (with power P <3 fW), we 
find that, at low temperature T, the resistance at the Dirac point Ro{H) undergoes a 1000-fold 
increase from ~10 kfl to 40 Mfi within a narrow interval of field. The abruptness of the increase 
suggests that a transition to an insulating, ordered state occurs at the critical field He- Results 
from 5 samples show that He depends systematically on the disorder, as measured by the offset gate 
voltage Vq. Samples with small Vo display a smaller critical field He- Empirically, the steep increase 
in Ro fits acccurately a Kosterlitz-Thouless-type correlation length over 3 decades. The curves of Ro 
vs. T at fixed H approach the thermal-activation form with a gap A ~15 K as // ^ H~ , consistent 
with a field-induced insulating state. 

PACS numbers; 73.63.-b,73.21.-b,73.43.-f 
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I. INTRODUCTION 

In graphene, the low energy states display a linear 
energy-momentum dispersion described by the Dirac 
Hamiltonian. The observation of the integer quantum 
Hall (QH) effect by Novoselov et al. [Ell,!! and Zhang 
et al. 0, B [1] has sparked intense interest in this novel 
2D (two-dimensional) system. In a strong magnetic field 
H , the states are quantized into Landau Levels. As a re- 
sult of the Dirac dispersion, the energy En of the Landau 
Level (LL) (of index n) increases with the flux density B 
as "v/S, viz. En = sgn(n) 2ehvpB\n\, where vf is the 
Fermi velocity, e the electron charge, and h is Planck's 
constant. The Hall conductivity is observed to be accu- 
rately quantized as axy = (4e^//i)[n -I- i] — ve^ /h, where 
e is the electron charge, 2-Kh is Planck's constant, and 
v the sublevel index. A key question is the nature of 
the ground state at the Dirac point. In intense i?, the- 
ory predicts interesting broken-symmetry states driven 
by exchange and interaction. These are characterized as 
quantum Hall Ferromagnetism (OTp_0Jl i, [M [HI, [ll , 
or excitonic condensation [l^, [Tj, [Tsrilq. These col- 
lective states imply the existence of field-induced phase 
transitions, but the experimental situation is rather un- 
settled. Moreover, the proposed [H, [13, dl existence of 
counter-propagating edge modes at the Dirac point has 
further enriched the theoretical debate. Is the high-field 
Dirac point a QH insulator or a QH metal? 

Recently, we reported [l^ that the resistance at the 
Dirac point i?o begins to increase steeply ai B — 10-12 T, 
suggesting a transition to an insulating state. However, 
the results left open several key questions. Because Rq 
increased by only 1-decade (to 0.2 MTl) ^1^, we could 
not establish that the high-field state is truly insulating. 
In graphene, the extreme sensitivity to thermal runaway 
has been highly problematical for researching its high- 
H ground state [1, [1, [l9|]. Adopting a low-dissipation 
technique to avoid self-heating, we have measured the 
divergence in Rq to 40 Mil (~ 1500 h/e^) in 3 samples. 



Remarkably, the divergence is accurately described over 
3 decades by the Kosterlitz-Thouless (KT) correlation 
length. The singular nature of the divergence provides 
strong evidence that a 2D transition to an insulating state 
occurs when B exceeds a critical field He- The systematic 
variation of He with | Vb | (the gate voltage needed to bring 
the chemical potential to the Dirac point) implies that 
disorder is very effective in delaying He to higher field 
values. In all samples investigated to date, the transition 
to the insulating state is reached in fields below 35 T. 



II. EXPERIMENTAL DETAILS 

Empirically, problems associated with self-heating in 
graphene become serious when the power dissipated P 
exceeds ~10 pW for bath temperature Tf, < 1 K. As dis- 
cussed in the Appendix, self-heating below 1 K leads to a 
number of spurious features caused by thermal instability 
in the sample. We adopted a simple voltage-controlled 
technique with ultra-low dissipation that avoids this dif- 
ficulty, and allows the divergence in Rq to be measured 
reliably to 40 Mil. An ac source maintains a fixed volt- 
age amplitude (40 /xV) across the sample in series with a 
lOO-kfi resistor (details are given in the Appendix). By 
phase-sensitive detection of both the current / and the 
voltage Vxx, we have made 4- probe measurements of i?o 
with ultra-low dissipation [P decreases from ~3 fW at 
10 T to 40 aW above 25 T). Moreover, for T <1.5 K, the 
sample is immersed in liquid "^He so that the electrons in 
graphene are in direct contact with the bath. The largest 
reliably-measured Rq is now 40 Mil (limited by the in- 
put impedance 100 Mil of the preamplifier). The new 
results provide an enlarged view of the interesting region 
in which Rq diverges. The samples K52 and J24 have off- 
set voltages Vo much larger than that in K7, the sample 
investigated in detail in Ref. [19]. All samples except J18 
were measured as-fabricated. Sample J18 was subject to 
a i hr. anneal in He gas at 80 C which decreased Vb to 
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8 V. However its large He suggests that its initial value 
of Vo (before annealing) is very large. 

In Samples K52, J18 and J24, the spacings a between 
voltage leads are 3.5, 2.75 and 3 /im, while the widths w 
are 3, 3 and 2 /im, respectively. 




the curves at 0.3 and 1.5 K (which cannot be distin- 
guished) show that i?o diverges to 40 MO, with a slope 
that steepens rapidly with H . The 3-decade increase (40 
kO 40 MO) occurs within the narrow interval 17-27 
T. We find that the observed divergence is too steep to 
fit a power-law of the kind i?o ^ {He — H)^", with a > 
and Hr a critical field. 
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FIG. 1: (color online) (Main Panel) Divergence of the resis- 
tance -Ro at the Dirac point with B at T = 0.3, 1.5, 5 and 
27 K (Sample K52). At 27 K, the increase in _Ro is quite 
moderate (to 190 kfi at H = 31 T). At T = 0.3 K, however, 
Ro exceeds 20 MQ above 27 T. The curves at 0.3 and 1.5 K 
undergo a 1000-fold increase (40 kfl to 40 MQ) in the narrow 
field interval 17-27 T. In high B, the 5 K curve deviates sig- 
nificantly from them. The inset shows the behavior of Ro vs. 
H in greatly expanded scale (xlOO). The voltage-regulated 
technique used for these measurements dissipates ~3 fW at 
10 T and 40 aW above 25 T. 



FIG. 2: (color online) Plot of Ro vs. l/\/l — h for Samples 
K52, J18 and J24 at T = 0.3 K, where h = H/Hc (the curve 
for K52 is displaced by 1 decade for clarity). Ro is expressed 
as sheet resistance (fi/D). For each sample. He is the opti- 
mal value that gives the best straight-line fit vs. l/Vl — h 
(curvature is noticeable if He is altered by ±0.1 T from this 
value). For Samples K52, J18 and J24, He equals 29.1, 32.1 
and 35.5 T, respectively. For K52, the thin solid line is the 
expression R^{h) = 440 exp[26/Vl - h], with b = 1.54. The 
data match R^{h) very well over nearly 3 decades in Ro. 



III. FIELD DEPENDENCE OF Ro 

Figure[T](main panel) shows curves of Rq vs. H in K52 
at temperatures T — 0.3 to 27 K. As evident in the curves 
at 0.3 and 1.5 K, i?o undergoes a very steep, divergent in- 
crease when H exceeds ~25 T. The region just before the 
divergence occurs is shown in greatly expanded scale in 
the inset. At 27 K, the increase in Rq is relatively mod- 
est (~20) between B — and 31 T (inset). However, 
as T decreases to 5 K, the increase steepens sharply, as 
reported [l^ for K7. Further cooling from 5 to 0.3 K 
changes the profile only very slightly. In the main panel, 



As in Ref. [l9|, we compare the divergence with that 
predicted for the Kosterlitz-Thouless (KT) transition. In 
2D systems described by the XY model, the ordered 
phase is destroyed at the KT transition by unbinding 
of pairs of topological excitations (e.g. vortices). As the 
transition here is induced by varying the applied field 
H, we replace the reduced temperature t by the reduced 
field h = H/Hc, with He the critical field. In the limit 
h ^ 1^ , the KT correlation length ^ diverges as 

f = aexp[6/Vl - h], (1) 

where a is the lattice parameter and b a number I. 

The quality of the fit to Eq. [T] is best revealed in a 
plot of logi?o vs. the quantity l/Vl — h. In Fig. [21 
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we have plotted logi?o in 3 samples K52, J19 and J24 
against l/\/l — h. In each sample, the value of He is 
adjusted to maximize the high-field portion of the plot 
that falls on a straight line (this is the only adjustment 
made). In Samples K52, J18 and J24, the inferred values 
of He are 29.1, 32.1 and 35.5 T, respectively. The values 
of Rq in K52 fit the straight line representing the expres- 
sion R^{h) = 440exp[2VVT^], with He = 29.1 T and 
b — 1.54. The value of b is consistent with the KT tran- 
sition. The 3-decade span is strong evidence that Eq. [1] 
accurately describes the divergence in Rq, and supports 
the inference that, at low T, we are observing a 2D KT- 
type phase transition to a high-field ordered state that is 
insulating. 
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FIG. 3: (color online) Comparison of the curves of _Ro vs. 
H at 0.3 K in log-log scale in 5 samples K7, K22, K52, J18 
and J24 with gate voltage offsets Vo = 1, -0.6, 3, 20, - and 
24 V, respectively (Vb before annealing is not known in J18). 
-Ro is expressed as sheet resistance H/D. In K7 and K22 
(which have small |Vo|), the divergence in Rq occurs at a 
lower H [l^- The femtowatt-dissipation technique applied to 
K52, J18 and J24 was not available for K7 and K22 (their 
curves were limited to _Ro <0.3 Mf2). The dashed curve is the 
fit of the K52 data to 7?^. 



we have obtained evidence that the zero-field mobility /ie 
is strongly suppressed if Vq is large. In Fig. [4l we display 
curves of R^x vs. the unshifted gate voltage in a batch 
of samples that includes K22. The curves are taken at 295 
K or 4 K (as indicated). In each sample, the width AVg 
of the peak at the Dirac point directly measures l//ie- 
Clearly, the width increases dramatically with Vq. These 
results support the inference that a large offset Vb gives 
rise to large electronic disorder which enhances disorder 
scattering and suppresses [i^- In turn, in an intense field, 
the transition field H^ is pushed to higher values. While 
these correlations do not exclude other factors that may 
influence H^, we have found that Vb is the single most 
reliable predictor of the field scale at which the divergence 
onsets at low T . 
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FIG. 4: Curves of the longitudinal resistance R^x vs. un- 
shifted gate voltage Vg in zero-H in a series of samples at 295 
K and 4 K (as indicated). The width of the peak in Rxx in- 
creases systematically with Vo (located by the peak position). 
This implies that the mobility fie is very low in samples with 
large Vo. 



The significant spread of He inferred from the fits in 
Fig. [5] is in accord with Ref. [l^ which reported that 
He correlates with the offset voltage Vq- In Fig. [2l the 
samples K52, J18 and J24 display a much larger critical 
field He than the sample K7 (with Vb = 1 V and He = 18 
T) studied in detail in Ref. [l^l- Figure [3] plots together, 
in log-log scale, the curves of Rq vs. H in the 5 samples 
for which we have detailed high-field transport results. 
The systematic shift to higher fields of the divergence 
(in the order K22, K7, K52 and J24) is matched by the 
increase in their Vq (-0.5, 1, 3, 20 and 24 V, respectively). 
The value of Vq before annealing in J18 is not known. The 
dependence of He on Vb is non-linear. Initially (for < 
Vb < 4V), He increases rapidly, but appears to increase 
rather slowly when Vq exceeds 20 V. 

From a study of how Vb affects the zero-field transport. 



IV. DOPING DEPENDENCE 

Further insight into the nature of the divergence is ob- 
tained by viewing the behavior of the longitudinal resis- 
tance Rxx vs. Vg in a narrow gate window around the 
Dirac point at fixed B (with T kept at 0.3 K). FigurelS^ 
displays a series of curves of R^x (in log scale) vs. Vg in 
K52 for fields 10 < < 31 T. At 10 T, Rxx displays 3 
well-separated peaks corresponding to n = LL at 20 V 
and the n = ±1 LL's at 6 and 38 V, respectively. For 
H >20T, the n = ±1 levels move out of the gate- voltage 
window. We focus on the divergent enhancement of the 
central peak as B increases to 31 T. The key feature 
is that Rxx at the Dirac point (Vg = 20 V) rises most 
rapidly especially for B >25 T. 
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FIG. 5: (color online) Variation of Rxx (Sample K52) vs. the 
gate voltage Vg in the interval (0< Vg < 38 V), with B fixed at 
selected values f0-3f T, and T = 0.3 K (Panel a). At B = f 
T, the central peak (n = LL) is well separated from the LL 
peaks labelled as n = ±1. At B = 20 T and higher, the n = 
±1 LL's fall outside the gate window. With increasing B, the 
central peak increases rapidly and broadens. At the maximum 
B (31 T), at the Dirac point {Vg = 20 V) diverges to 
values above 40 Mf2. The voltage-regulated technique has 
poor resolution when R^x falls below 0.3 kf2. Panel (b) shows 
the contour plot of Rxx{n2D, B) at 0.3 K in the n2D-B plane 
(color bar of Rxx shown on right). The dashed lines trace the 
sublevel degeneracy \/2-k(.%. The density is given by n2o — 
CVg/Ae where C and A are the capacitance and area of the 
device, respectively {C/A = 1.14 x 10"* Fm"^). 
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FIG. 6: (color online) The T dependence of Ro in Sample 
J18 with H fixed at selected values in the interavl 20< H < 
31.2 T. As H approaches He (32.1 T), Ro approaches the 
thermally activated form e^'''^ . The thin dashed curves are 
fits to this form with the gap A — 12.9 and 14.7 K at = 
30 and 31.2 T, respectively. At the highest field (31.2 T), Ro 
increases by a factor of 310 between 16 and 2 K. 



sensitive to slight deviations away from exact charge neu- 
trahty (Fig. [5]). As we decrease H below He, the ordered 
state is unstable to the spontaneous unbinding of (vortex- 
like) topological excitations which have a mean spacing of 
^ (Eq. [T]). Because Rq fits accurately to over 3 decades, 
we infer that the conductance scales as the density of ex- 
citations. Hence the excitations are charged, and they 
carry the entire current / in the limit T — > 0. As this 
conduction channel is qualitatively distinct from ther- 
mally activated carriers, it may account for the unusual 
"saturation" behavior of Rq vs. T [l^. As T ^ 0, with 
H fixed near He, the conduction crosses over at 1 K from 
a steep, thermally activated channel to a T-independent 
channel carried by the excitations (Fig. [Dd). 



It is instructive to view R^x (at 0.3 K) as a contour plot 
in the n2D~B plane where n2D is the 2D density of carri- 
ers doped by gating (Fig. [SJd). The color bar (right) gives 
the magnitude of Rxx- Interestingly, the steep increase in 
Rxx with B is confined to the region between the dashed 
lines, which trace the sublevel degeneracy l/27r^^. This 
suggests that only the states within the lowest sublevels 
(on either side of = 0) are affected by the opening of the 
gap A. The contours appear to converge to a rounded 
cusp at n2D = 0, but with a curvature that increases 
rapidly with B. At the largest Rxx (white region), the 
contour resembles a a narrow, sharp wedge. The contour 
pattern suggests that He increases very rapidly from the 
value 29.1 T, as \n2D\ deviates from 0. 

The physical picture implied by the results is that, at 
the Dirac point, a field-induced transition to a gapped, 
insulating state occurs at He- The value of He is highly 



V. TEMPERATURE DEPENDENCE 

The T dependence of i?o in J18 with H fixed at selected 
values is displayed in semilog scale in Fig. [S) At the rel- 
atively high field H = 20 T, the 2-fold increase in Rq 
between 16 and 2 K is quite modest. However, above 20 
T, the T dependence steepens rapidly. As H is increased 
towards the critical He (=32.1 T), the profile of Rq moves 
ever closer to the activated form Rq ~ e^/'^ . We have 
plotted fits to the activated form (dashed curves) at the 
2 highest fields = 30 and 31.2 T. At the highest H, the 
measured i?o tracks closely the dashed curve until T falls 
below 3 K where it deviates downwards. As H decreases 
further from He, the deviations start at a higher T. The 
gap value A equals 12.9 and 14.7 K at 30 and 31.2 T, re- 
spectively. The latter provides an estimate of the energy 
gap in the ordered state for J18. Establishing that Rq is 
thermally activated when H > He is an essential goal, as 
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it shows that the ground state above He is a true insula- 
tor with a well-defined gap order parameter (as opposed 
to a state in which the carriers are strongly localized). 
Although the curves in Fig. [Hlcome close to establishing 
this result, measurements of Rq vs. T in samples with 
more accessible He are desirable. 

Interestingly, throughout the pre-transition region 
(20 < H < He), Rq also displays large increases with 
decreasing T. The carriers are strongly affected by the 
impending insulating state. We interpret these changes 
as reflecting very strong fluctuations in the order pa- 
rameter that characterizes the ordered insulating phase. 
The strong T dependence is also apparent in the fixed-T 
curves shown in Fig. [T] 

Quite apart from the thermal activation argument, 
there are other evidence to suggest that the divergence 
is not consistent with electron localization. As evident 
from measuring the widths of the peaks of R^^. vs. Vg 
taken in zero B, samples with smallest |Vo| are the least 
disordered. The electron mobility fj,e decreases from ~2.5 
to 0.5 T-i as IVol increases from 0.5 to 20 V. With this 
trend in mind, we compare in Fig . [3] the profiles of Ro 
vs. B in K52 with K7 and K22 [ig. In Samples K22 and 
K7 (with Vq = -0.6 and 1 V, respectively), the divergence 
in Rq is apparent in relatively low B (below 12 T). By 
contrast, we must go to much higher fields (>20 T) in 
K52. The dashed line is the fit to R^ in K52 {H > 18 T) 
described above. 

In the localization scenario, the observed divergence 
of Rq in strong field is explained by postulating that B 
induces localization of the electrons. However, applying 
this reasoning to the 3 samples in Fig. [31 we would con- 
clude that a modest B is sufficient to trigger the localiza- 
tion in clean samples, but very intense fields are needed 
in dirtier samples. This implies that disorder and field 
act in opposition to bring about localization, which is in 
conflict with physical intuition. In addition, localization 
induced by B cannot lead to the singular divergence ob- 
served in Ro (Fig. For these reasons, we believe that 
localization is not a viable explanation for the divergence 
in Rq. 



VI. DISCUSSION 

In the absence of Zeeman splitting and electron- 
electron interaction, the n = LL is 4-fold degenerate 
corresponding to the (physical) spin degeneracy indexed 
by a = ±1 and the K and K' valley degeneracy indexed 
by r = ±1. At v = 0, the energy of the n — LL, 
is zero. The effect of interaction in producing a broken- 
symmetry ground state has been investigated by several 
groups. To discuss our experiment, it is convenient to 
distinguish 2 different theoretical scenarios for the n = 
LL. 

In one scenario, the Quantum Hall Ferromagnet (QHF) 
models, the exchange energy E^x ~ VB leads to ferro- 
magnetic polarization of the physical spins 0, EMM 



llTj l . This produces a spin gap in the bulk without affect- 
ing the valley degeneracy, i.e. E"^ = a{iiBB + E^x) with 
jiB the Bohr magneton. Near the edge of the sample, the 
residual valley degeneracy is lifted by the edge potential. 
An important consequence of the QHF scenario at = 
is the existence of spin-filtered counter-propagating edge 
(CPE) modes which result in a residual conductance of 
2e^//i regardless of the magnitude of the spin gap in the 
bulk [12 [13, [13 ■ [In principle, the CPE modes are not 
present if the exchange polarizes instead the valleys to 
produce the bulk gap (this involves the same exchange 
energy E^x)- However, most investigators favor the spin- 
polarization scenario in graphene because it is augmented 
by the Zeeman energy fj-sB (the valley-polarization sce- 
nario is also called the QHF following the original us- 
age [20|] in quadratic, bilayer GaAs-based devices).] 

In the second scenario, called magnetic catalysis [ll|, 
[Tsl . [3, [H, [3|, the field component B± normal to the 
graphene sheet triggers electron-hole condensation. The 
instability introduces a mass term to the Dirac equa- 
tion which leads to the order parameter [l^, [3, Il5l | 

(^*) = J^rail'^rAal'^ " IV'tB^P), whcrC VrAa and ^JrBa 

are the wave functions of electrons of spin a and val- 
ley T at sites A and B, respectively. The instability - 
a solid-state realization of chiral-symmetry breaking in 
(2+l)D [2l[ - results in preferential occupation of, say, 
the A sublattice sites over the B sites, and drives the sys- 
tem into an insulating state. Significantly, the instability 
is strongest for n = 0. 

The steep increase in i?o vs. H first reported in Ref. 
[Tot implies that at large H, the ground state at the Dirac 
point has a resistance at least 20 times larger than the 
quantum hje^. Although the measurements were limited 
to i?o < 0.3 Mil, the upturn appeared to diverge at a 
critical field H^ suggestive of a singular field dependence. 
The findings are clearly at odds with the existence of 
CPE modes (see, however, the results in Ref. [IS]). 

In the present report, we have extended by a factor 
of 200 the range of resistance measurements and shown 
that, at 0.3 K, the increase in Rq is truly divergent as well 
as singular. Moreover, this behavior has been observed 
in all samples investigated to date by us in high fields. 
The evidence amassed clearly establish that the high- 
field ground state at the Dirac point is a true insulator 
(at least for samples prepared on a Si02 substrate). Ipso 
facto, the CPE modes do not exist in the insulating state. 
However, our results do not preclude them at low fields. 

Lately, several groups have considered how the CPE 
modes are affected by intense field. A very interest- 
ing possibility is that an intense field destroys the CPE 
modes in a field-induced transition. It has been pointed 
out to us that the CPE modes are not protected against 
2-particle exchange scattering with spin flip [2^. As the 
exchange energy increases with _B, the increased scatter- 
ing rate could lead to a gap in the edge modes. 

In the magnetic catalysis scenario, Gorbar et al. [lB| re- 
cently considered the competition between the mass gap 
(^'5') and the spin gap (augmented by Zeeman energy) 
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and inferred that CPE states exist only above a critical 
field Bcr- 

Shimshoni et al. [23] have proposed that scattering off 
magnetic impurities can lead to strong localization of 
electrons in CPE modes which can mimic a KT transi- 
tion. However, this scenario needs to be reconciled with 
the observed nearly activated behavior of Rq as well as 
the variation of with Vq in different samples. 



VII. APPENDIX 

Near the Dirac point, resistance traces are strongly 
distorted when the Ohmic heating P exceeds 10 pW at 
bath temperatures below 1 K. As examples, we plot 
in Fig. [7| resistance traces (with / fixed). In Panel a, 
the inferred curve of Rxx vs. Vg, measured with B — Zl 
T, / = 10 nA and T — 0.3 K, shows a pronounced dip 
near the Dirac point caused by self-heating instability 
(the true i?o exceeds 10 Mfi). Panel b shows "i?o" vs 
B measured at fixed /. In the curve for K23 (at 5 K), 
self- heating reverses the trend of Rq. The downturn is 
avoided when / is decreased to 1 nA, but the measured 
curve (in K22 at 0.3 K) is still greatly suppressed from 
the true divergent profile. 
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FIG. 7: (color online) Spurious features caused by sample 
self-heating in graphene. Panel (a) shows a gate-sweep mea- 
surement at T = 0.3 K of Rxx vs. Vg in K52, with / fixed 
at 10 nA (dc). At B = 31 T (red curve), severe self-heating 
inverts the peak at the Dirac point (as shown in Fig. [T] Ro 
actually exceeds 40 Mfi). Heating effects are less severe in 
the 20-T curve. Panel (b) displays curves of "Ro" vs. B mea- 
sured in Samples K22 (at 0.3 K) and K23 (at 5 K) with I 
fixed at 1 and 10 nA (dc), respectively. When Ro > 0.2 Mfl, 
self-heating produces the spurious shoulders and broad peaks, 
whose positions and shapes depend on I. 



Figure [5] is a schematic of the measuring circuit em- 
ployed in the ultralow-dissipation technique. A nomi- 
nally constant ac voltage (~40 //V) of frequency 3 Hz is 
applied across the sample in series with a 100-kJ7 buffer 
resistor. The current passing through the sample is mea- 
sured by a Keithley picoammeter whose output is phase- 




FIG. 8: (color online) Schematic of the low-dissipation, 
voltage-regulated circuit used in the experiment. Lock-in 
(amplifier) A produces a regulated voltage emf (3 Hz) that is 
reduced to an amplitude of 40 /iV by a 100:1 voltage divider. 
The signal goes through a tt filter, a buffer resistor (100 kf2) 
and a low-pass filter before entering the dewar. The AC cur- 
rent passing through the graphene sample is measured by the 
picoammeter (Keithley), whose output is phase-detected by 
Lock-in A. The longitudinal voltage Vxx and Hall voltage Vxy 
are phase-detected by Lock-ins Vi and V2, respectively, after 
transmission through a bank of 7r-filters and high-impedance 
(100 MSI) pre-amplifiers. As shown, all wires entering the de- 
war are buffered by 1 kSl nichrome thin-film ceramic resistors. 
The inset (lower left) shows Sample J 18. 



detected by the lock-in amplifier A. Simultaneously, the 
longitudinal voltage Vxx and Hall voltage Vxy are phase- 
detected by 2 other lock-ins. As shown, all wires entering 
the dewar are filtered and buffered to exclude extraneous 
RF signals which may be a potential source of sample 
heating. 

Using the ultralow-dissipation technique, we com- 
pletely avoid the thermal runaway problems illustrated 
in Fig. [T] At selected fields, we have performed I-V 
measurements to check that self-heating is not skewing 
the results even at our lowest T (0.3 K). Figure O shows 
curves of / vs. Vxx at the Dirac point in K52 at T = 0.3 K 
with H fixed at 20 and 24 T. The linearity implies that 
self-heating is not observable up to a bias voltage of 2 
mV. Since all the curves displayed in the main text were 
taken with a bias of 40 /iV, we are comfortably within 
the Ohmic regime. 

We thank H. A. Fertig, P. A. Lee, D. Abanin, V. 
A. Miransky, E. Shimshoni, F. D. M. Haldane and D. 
N. Sheng for valuable comments, and acknowledge sup- 
port from NSF-MRSEC under Grant DMR-08 19860, and 
from the Princeton Center for Complex Materials. The 
experiments were performed at the National High Mag- 
netic Field Laboratory, which is supported by NSF Co- 
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FIG. 9; (color online) The current-voltage curves measured 
at the Dirac Point {Vg = Vo) in Sample K52 at T = 0.3 K at 2 
the fields H ^ 2Q and 24 T (using the circuit in Fig. [H|. The 
linearity of I vs. Vxx implies that, at these fields, self-heating 
is not observable at these power-dissipation levels (the dashed 
lines are guides to the eye). 
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